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RESUMEN 



El resumen sera traducido al espanol por los editores. Starburst galaxies are generally associated with 
extended X-ray and radio halos, giving a clear hint of an outflow of gas and relativistic particles from the disk 
fT^ , into the halo. The driving agents are, not surprisingly, active star forming regions, injecting hot gas and cosmic 
' rays generated by supernova remnant and superbubble shock waves. The dynamics of the outflow and the 
thermal evolution of the plasma are strongly coupled, and therefore a self-consistent model is necessary for a 
satisfactory description. Since the plasma temperature is between one and a few million Kelvin, the halo is most 



^ . conspicuous in soft X-rays. It will be shown that X-ray data obtained recently with XMM-Newton EPIC pn 
bear clear spectral signatures that are in strong disagreement with an isothermal halo in coUisional ionization 
equilibrium (CIE). Instead a temperature structure in X-ray halos is observed. It will be demonstrated that 
' a galactic wind outflow model in which the non-equilibrium ionization structure is calculated self-consistently 

can give a satisfactory and physical explanation for the X-ray spectral characteristics, e.g. in case of the local 
starburst galaxy NGC 3079. In particular, high subsolar abundances, that have been reported from X-ray 
observations in some cases, are shown to be artefacts of CIE spectral fit models. It is found that spectral 
' models are strongly constrained by the presence of diagnostic lines such as O VII and O VIII as well as Fe-L 
. line complexes. Thus a detailed spectral X-ray model will help to provide important galactic wind parameters, 
' such as mass loss rate and velocity profile, which may also serve as valuable input for galactic outflows in the 
early universe. 



ABSTRACT 

Starburst galaxies are generally associated with extended X-ray and radio halos, giving a clear hint of an outflow 
^ ^, of gas and relativistic particles from the disk into the halo. The driving agents are, not surprisingly, active 
star forming regions, injecting hot gas and cosmic rays generated by supernova remnant and superbubble shock 
waves. The dynamics of the outflow and the thermal evolution of the plasma are strongly coupled, and therefore 
a self-consistent model is necessary for a satisfactory description. Since the plasma temperature is between 
^3 . one and a few million Kelvin, the halo is most conspicuous in soft X-rays. It will be shown that X-ray data 
obtained recently with XMM-Newton EPIC pn bear clear spectral signatures that are in strong disagreement 
with an isothermal halo in coUisional ionization equilibrium (CIE). Instead a temperature structure in X-ray 
' halos is observed. It will be demonstrated that a galactic wind outflow model in which the non-equilibrium 
ionization structure is calculated self-consistently can give a satisfactory and physical explanation for the X- 
ray spectral characteristics, e.g. in case of the local starburst galaxy NGC 3079. In particular, high subsolar 
abundances, that have been reported from X-ray observations in some cases, are shown to be artefacts of CIE 
spectral fit models. It is found that spectral models are strongly constrained by the presence of diagnostic lines 
such as O VII and O VIII as well as Fe-L line complexes. Thus a detailed spectral X-ray model will help to 
provide important galactic wind parameters, such as mass loss rate and velocity profile, which may also serve 
as valuable input for galactic outflows in the early universe. 
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1. INTRODUCTION in a state of enhanced star formation (Searle et al. 

r,^, , ,. r 11 1 i J 1 i i 1973). These so-called starbursts are intermittent 
Ihe observation oi optically selected late-type ^ . ' , , . . , , . r • ,, 
, . 1 J i J.V, 1 • j-u J- 4.1, ui 4. u during galaxy evolution with a duration of typically 
galaxies led to the conclusion that the bluest ob- « , , , , , , r 
. , u 4. 11 I, 1, 4. 1-11 10 yrs. It has been argued that the total rate of 
jects could not all be young, but more likely are , . , r . . , . , , , . , 
high-mass star formation withm these short (with re- 

'^Max-Planck-Institut fiir extraterrestrische Physik, Garch- Spect tO a Hubble time) periods is comparable in the 

ing, Germany. 
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local universe to the rate in spirals during the whole 
time of quiescent star formation (Heckman 1997). 
Undoubtedly, the effects of a starburst on the ISM 
are dramatic. A high supernova rate disturbs and 
eventually disrupts the gaseous disk and supernova 
remnant (SNR) heated gas and accelerated energetic 
particles (cosmic rays) are injected into the base of 
the galactic halo. 

It has been shown theoretically, that the com- 
bined pressure gradients of gas, cosmic rays (CRs) 
and MHD waves are able to drive a galactic wind 
even in cases of moderate star formation, like in our 
own Galaxy (Breitschwerdt et al. 1991). A nec- 
essary condition is that there is a dynamical cou- 
pling between the plasma and the CRs. Such a cou- 
pling can be provided by strong resonant scattering 
off self-excited MHD waves, which are generated by 
a large scale CR gradient pointing away from the 
galaxy and thus induce a so-called streaming insta- 
bility (e.g. Kulsrud & Pearcc 1969). On the other 
hand, the overpressure in supcrbubbles generated by 
rich OB clusters is sufficient to drive a largely ther- 
mal outflow. In case of starburst galaxies the con- 
ditions are even more extreme, and in galaxies like 
M82, the galactic wind (or superwind, as it is some- 
times called) is purely thermally driven (Chevalier & 
Clegg, 1985; Volk et al. 1996). 

In recent years the study of starburst phenomena 
and galactic winds has received a tremendous boost. 
This was largely driven by observations, notably the 
availability of 8 m-class telescopes such as the VLT 
and Keck with superb instrumentation, as well as 
the HST. Selecting galaxy samples by colour and 
measurement of the UV flux of galaxies indicated a 
peak in the star formation rate (SFR) near redshift 
z ~ 1.5 (Madau et al. 1996), while in more recent 
observations it was argued for a near constancy be- 
tween z ~ 1.5 — 4 (Steidel et al. 1999). A major 
uncertainty in these studies was the amount of ex- 
tinction in the UV. Therefore both extinction and 
redshift should be included in the SFR analysis. A 
detailed study of the northern Hubble Deep Field 
with HST NICMOS, in which both the SFR was de- 
termined individually for galaxies from the 1500A 
UV flux and a spectral energy distribution template 
fitting method was applied to determine redshift as 
well as extinction, yielded the result (Thompson ct 
al. 2001) of an increase in SFR between z = 1 — 2, 
then a fall-off from 2 = 2 — 3 and a plateau region be- 
tween z = 3 — 5 (see Fig.^; as a clear result, the SFR 
was higher up to redshifts of ~ 5 than at present. 
Naturally, an enhanced SFR implies an enhanced 
starburst activity during the infancy of galaxy evo- 
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Fig. 1. Star formation rate as a function of redshift in 
a comoving volume including also extinction (Thompson 
et al. 2001) 
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Fig. 2. Lya emission line from the outflow of the star- 
forming galaxy J123649.2-F621539 (Dawson et al. 2002). 
The fitted line profile could be modeled by 3 components: 
a large amplitude narrow Gaussian for the recombination 
of H"*" , a small amplitude broad Gaussian (redshifted by 
320 km/s) for backscattered red-wing emission of far-side 
Lya photons of the expanding wind, and a broad Voigt 
absorption profile (blueshifted by 360 km/s) to model the 
blue decrement due to near-side absorption of photons. 

lution and galactic winds are expected to have been 
fairly common at that time. Such an assumption is 
indeed confirmed by observation. A recent serendip- 
itous discovery of a galactic wind at z = 5.19 by 
Dawson et al. (2002) revealed a strong asymmetric 
Lya emission line (sec Fig. The best fit model 
line profile resulted in a red- and blueshifted wing, 
indicating an outflow velocity of 320 — 360 km/s. 
As will be shown below, such velocities are entirely 
consistent with wind speeds derived from modeling 
the X-ray emission of local starburst galaxies. An 
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important consequence of winds is the pollution of 
the intergalactic medium (IGM) with chemically en- 
riched material. In addition to metals, also entropy 
is added to the IGM. 

This may be an explanation for the so-called 
"entropy floor" postulated by comparing the X-ray 
luminosity versus temperature relationship. Going 
from massive clusters of galaxies to poor groups there 
is a deviation from a power law towards the low mass 
end. This has been interpreted in terms of an "en- 
tropy floor" that dominates the gravitational heat- 
ing of the poor clusters and groups as the potential 
wells become shallower (Ponman et al. 1999). Thus 
the preheating of the IGM by starburst driven galac- 
tic winds starts to dominate the release of gravita- 
tional energy. Also, the fairly high abundances of 
Z - 0.3 Zq (e.g., Molendi et al. 1999) found in 
the intra-cluster gas may be explained by ejection of 
metals from galaxies by massive winds during initial 
starburst phases. 

In the local universe, evidence for galactic winds 
mainly stems from observations of starburst galax- 
ies, such as M 82 or NGC 253. In these more evolved 
objects, starbursts are thought to be triggered by a 
substantial disturbance of the gravitational poten- 
tial, e.g. by interaction with a companion. Never- 
theless, as it has been pointed out, local starbursts 
are similar in origin, evolution, duration and its ef- 
fects on the IGM to those at high redshifts (e.g. Pet- 
tini et al.2000). If one accepts these findings as a 
working hypothesis, nearby starburst galaxies would 
serve as ideal laboratories for studying the effects of 
starbursts on the early evolution of galaxies and the 
interaction with the IGM. Since the starburst driven 
outflows are mostly thermal, they are best observed 
in soft X-rays. 

2. OBSERVATIONS OF LOCAL STARBURST 
GALAXIES 

Direct evidence for supernova (SN) heated galac- 
tic halos stems from imaging and spectroscopy of dif- 
fuse soft X-ray emission. This has been impressively 
demonstrated by a recent XMM-Newton observation 
of NGC 253 (Pietsch et al. 2001). From Fig. it 
can be clearly inferred that (i) the diffuse emission 
is soft (cf. 0.2 — 0.5 keV band colour coded in red), 
and (ii) that the outflow is extended and distributed 
over the whole disk rather than being constrained 
solely to the nuclear region. The nearby (distance 17 
kpc) edge-on (inclination to the line of sight about 
85°) spiral galaxy NGC3079 has been observed with 
XMM-Newton for a total of 25 ksec. The galaxy is 
classified as an SBc LINER with distinct nuclear ac- 
tivity. The object is an ideal target for analyzing the 




Fig. 3. The extended soft X-ray emission of the nearby 
edge-on starburst galaxy NGC 253, as observed by the 
EPIC pn and MOS cameras onboard XMM-Newton in a 
three-colour image (red: 0.2 - 0.5 keV, green; 0.5 - 0.9 
keV, blue: 0.9 - 2.0 keV; overlayed contours represent the 
2.0 - 10.0 keV band; Pietsch et al. 2001). For reference 
the D25 ellipse is plotted to show the optical extension of 
the gas disk. The inlay is a zoom-in of the nuclear region, 
with ROSAT detected sources marked by squares. 

diffuse X-ray emission in the disk and the surround- 
ing galactic halo, since galactic foreground absorp- 
tion is low (column density Nh < lO^^cm"^). The 
target was chosen to investigate both the problem of 
starburst AGN connection, and to see if there was a 
wide-spread extranuclear X-ray emission, for which 
we wanted to derive its spectral properties. The ob- 
servation was carried out as part of the XMM SSC 
Guaranteed Time Program. Rejection of high back- 
ground times led to an overall useable integration 
time of 16.4 ksec. We detected a huge soft (0.2 — 1.0 
kcV) X-ray halo (s. Fig. ^ extending perpendicular 
to the galaxy disk to about 17.5 kpc (Brcitschwcrdt 
et al., 2003). The emission is therefore expected to 
be dominated by thermal lines, which are abundant 
in the soft band. The spectrum of a conspicuous 
emission region indeed proves this point. Moreover 
it clearly shows that the galactic starburst must drive 
a thermal outflow, since we have strong indications 
for collisionally excited oxygen and iron L line com- 
plexes in the spectrum (s. Fig. (T)). A further ar- 
gument for a halo connected to starforming regions 
in the underlying gaseous disk rather than to AGN 
activity is the softness of the halo, which completely 
disappears in the higher energy band (1.0 - 2.0 keV). 
The morphology of the halo also shows extranuclear 
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Fig. 4. The extended soft X-ray halo of the nearby edge- 
on starburst galaxy NGC3079, as observed by the EPIC 
pn camera onboard XMM-Newton in the 0.2 - 1.0 keV 
energy band (Breitschwerdt et al. 2003). Point sources 
are shown by circles. For reference the D25 ellipse is 
indicated to show the optical extension of the underlying 
gas disk. 

spurs supporting the starburst connection. 

3. SELF-CONSISTENT DYNAMICAL AND 
THERMAL MODELLING OF OUTFLOWS 

The injection of SN heated gas into the base of 
the halo generates a huge pressure gradient with re- 
spect to the ambient IGM. Consequently, the hot 
plasma can escape the gravitational potential well 
and expand to infinity, provided that the external 
pressure is sufficiently low. Otherwise, the expand- 
ing hot gas will become thermally unstable and rain 
down in the form of intermediate velocity clouds 
on the galactic disk (galactic fountain). This is 
to be expected if the radiative cooling time scale, 
Tcooi — 3nA{T) / {ksT) , is much less than the dy- 
namical flow time scale, Tf{z) = f^^dz'/u{z')] in a 
steady state flow this condition has to hold at least 
at the sonic point. Further out the flow is super- 
sonic and cannot be influenced by the source re- 
gion; mass loss is then inevitable. In a starburst 
driven flow this is certainly fulfilled as our calcu- 
lations show. The model as applied to an edge-on 
galaxy is sketched in Fig. El We have tested the 
hypothesis of a thermally driven superwind by mod- 
elling the outflow with our galactic wind hydrocode, 
in which the radiative cooling and X-ray emission 



Outflow (galactic fountain, wind, ...) 




Edge-on Galaxy 

Fig. 5. Sketch of the galactic outflow model. The wind 
is mainly driven by thermal pressure, but cosmic rays 
and MHD waves may also assist and are included self- 
consistently. 



is treated self-consistently with the dynamics in full 
non-equilibrium (NEI). In essence this means that 
we follow the time-dependent ionization structure of 
the outflow instead of using the usual assumption of 
collisional ionization equilibrium (CIE). The neces- 
sity in doing so has been been extensively discussed 
in Breitschwerdt & Schmutzler (1999). For brevity 
we just repeat here the main arguments. Firstly, 
the assumption of CIE is strictly never fulflUcd. A 
detailed balancing of the collisional ionization rate 
by the recombination rate would require three body 
collisions (one electron recombining with an ion, an- 
other one carrying away the excess energy), which 
are extremely rare in an optically thin and dilute 
plasma. Therefore collisional de-excitation is negli- 
gible and the coronal approximation holds. Instead 
radiative cooling, by which the excess energy is re- 
moved by a photon, is dominating by far. Except 
for resonant processes, the absorption probability of 
these photons is very low and they escape. But 
this means that both collisional ionization and ra- 
diative recombination are cooling processes for the 
plasma, which will hence be driven out of equilib- 
rium (Shapiro & Moore, 1976; Schmutzler & Tschar- 
nuter, 1993). CIE can therefore only hold as an ap- 
proximation during a limited period or due to some 
form of heat input. Secondly, it has been shown 
that if the plasma is in a dynamical state (e.g. ex- 
pansion), the dynamical time scale can be shorter 
than any atomic time scale and the plasma is char- 
acterized by a state of delayed recombination (Bre- 
itschwerdt & Schmutzler, 1994). Thirdly, there ex- 
ists a feedback mechanism between thermal and dy- 
namical evolution. The latter changes the density. 
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pressure etc. of the flow, which in turn changes the 
time-dependent ionization structure and thereby the 
radiative energy losses, which crucially determine the 
dynamics again. An approximate and efficient way 
in solving the hydrodynamics self- consistently along 
with the non-equilibrium ionization (NEI) structure 
is to calculate the temporal evolution separately 
and couple them by a time- dependent cooling rate 
C[n{t),T(t); Z], where n, T and Z are the density, 
temperature and chemical abundances, respectively 
(Breitschwerdt & Schmutzler, 1999). This iterative 
procedure converges fairly rapidly. The NEI code 
treats the time-dependent atomic processes (colli- 
sional ionization, recombination, autoionization, di- 
electronic recombination, to name just a few; for de- 
tails see Schmutzler & Tscharnutcr, 1993) of the ten 
most abundant elements and their respective ioniza- 
tion stages, and includes at present 1155 lines. The 
hydrocode has been described in detail by Breitschw- 
erdt et al. (1991). It entails a realistic gravitational 
potential of the galaxy, and the flow is calculated 
in steady-state in a flux tube geometry, which is a 
fair representation of the transition from plane par- 
allel flow close to the disk to spherical divergence 
far away. We expect that steady-state is justified 
during the starburst and as long as there are no 
dramatic changes at the inner boundary conditions. 
Fully time-dependent simulations, using the STAR- 
BURST99 code, are underway (Dorfi & Breitschw- 
erdt 2003), and conffim this by and large, but also 
show that abrupt changes in the SN rate lead to 
propagation of shocks, reheating of the halo and par- 
ticle acceleration out there to energies well in excess 
of 10^^ eV reached in individual SNRs. This has to 
be borne in mind when interpreting the spectra. 

In order to compare the simulated spectra di- 
rectly with the observations, we have folded them 
through the EPIC pn instrumental response. Tak- 
ing the instrumental field of view of ~ 30' x 30', 
we have projected the outflow cone at the galaxy's 
distance (see Fig. El onto it. Thus we do not only 
reproduce the form of the spectrum, but also the 
absolute count rate. We even went one step further, 
and instead of comparing the observed and synthetic 
spectra by eye, we have used the latter as input for 
XSPEC to actually fit the observed spectrum. The 
result shows (s. Fig.El) that the halo indeed exhibits a 
" multi-temperature structure" , fully consistent with 
an NEI outflow, in which supernova heated gas is 
injected into the halo at temperatures of 3.6 x 10^ 
K and density 5 x 10~^ cm~^ at an initial velocity of 
about 220 km/s (see Fig. [SI to give a satisfactory fit 
with Xred ~ A closer inspection of Fig.|7|reveals 
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Fig. 6. Dynamically and thermally self-consistent out- 
flow model of the starburst driven galactic wind of 
NGC3079 (Breitschwerdt et al. 2003). The solid line 
is the outflow velocity u{z), va{z) is the Alfven veloc- 
ity (dash-dotted line), c*{z) is the generalized speed of 
sound (dashed line) and VeBc{z) (dotted line) is the escape 
speed from the gravitational potential of the galaxy. 
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Fig. 7. Fit (x^ = 1.2) of the Epic pn observed spectrum 
(crosses = error bars) of the soft X-ray halo emission 
of NGC3079 between 0.2 - 2 keV obtained with a self- 
consistent non-equilibrium emission outflow model (solid 
line; Breitschwerdt et al. 2003). 

two distinct humps at 0.6 keV and 0.9 keV, caused 
by the O VII, O VIII and Fe-L line complexes, re- 
spectively. We emphasize that a single temperature 
CIE model, as e.g. by Raymond & Smith (1977) or 
a so-called MEKAL model are not able to reproduce 
the spectrum in a statistically acceptable way. The 
reason is, because CIE models are extremely temper- 
ature sensitive, each hump in our spectrum would 
require a different CIE temperature. In an attempt 
to achieve a fit, although not a very good one, a com- 
mon error is to decrease the element abundances to 
extremely low values (e.g. a few percent solar, as 
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been advocated in some publications) in order to get 
rid of prominent lines which would enforce a specific 
excitation temperature. We consider this as a pure 
artefact of the fit procedure, and moreover as un- 
physical, because an outflow that is mass loaded by 
chemically enriched material from the starburst re- 
gion should have rather high than low metallicities. 

4. CONCLUSIONS 

Star formation is a highly nonlinear process, and 
it has been known for some time, that for exam- 
ple it can be triggered by nearby SN explosions. 
On the other hand the rate cannot increase indef- 
initely, because the exhaustion of gas fuel renders 
it more and more inefficient, thereby forcing it into 
a self-rcgulating cycle. The highest rates known to 
date are realized by galaxy-galaxy interactions (with 
mergers as an extreme form), and starbursts as a 
common manifestation. In particular in the early 
universe when the average distance between galaxies 
was much smaller, starbursts could dominate their 
integrated luminosity. Both local and distant star- 
burst galaxies show that the enhanced star forma- 
tion rate is accompanied by a thermally driven out- 
flow, which can be traced to large distances in the 
soft X-rays. Although only a few percent of the en- 
ergy released is radiated in this wavelength range, 
the presence of lines of highly ionized an abundant 
species such as oxygen and iron in the spectrum 
serves as a clear fingerprint for the generation of a 
hot plasma and its dynamical evolution as a galac- 
tic wind. Therefore spectral modelling of outflows 
provides important quantitative information of mass 
loss rates, chemical enrichment of the IGM, and, pos- 
sibly a contribution to the WHIM (warm hot inter- 
galactic medium) that has been invoked to explain 
the so-called missing baryon problem (Cen & Os- 
triker 1999). Last but not least, the effect of a star- 
burst induced superwind on the further evolution of 
the host galaxy should not be underrated. Venting 
hot material away removes a lot of entropy thus en- 
suring a continuation of star formation in a more 
quiescent fashion in the disk. 
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